Lattice model calculations on aqueous acetone and tetramethyl urea 
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Self-organization of aqueous acetone and tetramethylurea, caused by lyotropic phase transfor- 
mations, is studied in frames of the Ising model. Using lattice Monte Carlo simulations, the con- 
centration dependent mixing schemes in these solutions are confirmed and the molar fractions cor- 
responding to the transitions between the mesophases of these mixing schemes are theoretically 
calculated. 

PACS numbers: 



INTRODUCTION 

In our study, we select two mixtures, aqueous solutions 
of acetone and tetramethyl urea, where hydrogen bond- 
ing or the longer range dipole-dipole coupling plays an 
illustrative role in revealing structure properties of wa- 
ter. 

Both these co-solvents are miscible with water over 
the entire concentration range and at all temperatures 
between RT (R is the gas constant) and boiling points. 
There is no lyotropic micellar phases, as the co-solvent 
molecules are rather small and globular in shape. Still, 
according to thermodynamic data, three composition re- 
gions can be identified in these mixtures, in each of which 
the mixing scheme is qualitatively different from those in 
other regions. The structures and interactions within 
these regions are certainly defined by the features of the 
constituent molecules, and can be roughly expressed in 
the following way. 

In the region of low concentration of the co-solvent, 
a small number of AC molecules are bound together by 
C=0 dipole attraction more tightly than in pure acetone, 
and for TMU solution, molecules are bound together by 
H = CHz dipole attraction, respectively. In solution, 
such a cluster enhances the hydrogen bond network of 
H2O in its immediate vicinity, with concomitant reduc- 
tion in the hydrogen bond probability of bulk H2O away 
from solute. So, such a regions were introduced as con- 
cerning to Mixing Scheme I [lj]. In the most surfactant- 
rich region, AC/TMU molecules form clusters in ex- 
actly the same manner as in pure liquid AC/TMU, in 
which AC/TMU molecules are bound together not only 
by C = O/H = CH3 dipole but also van der Waals inter- 
action. To the clusters of pure AC/TMU, H2O molecules 
interact as a single molecule. Hence, there is no hydrogen 
bond network of liquid H2O. This is Mixing Scheme III 
l] . Mixing Scheme II is operative, in which the solution 
consists of two kinds of clusters each reminiscent to I and 
III. The transition from Mixing Scheme I to II is associ- 
ated with an anomaly in a thermodynamic quantity, and 
it occurs in a narrow mole fraction range. 



For TMU, the effect of H 2 was studied 0. TMU 
molecules directly participate in forming the hydrogen 
bond network of H2O, presumably via H donor on NH2 
and H acceptor on C — O, up to almost the saturation. 
Hence, the hydrogen bond connectivity is not altered by 
the presence of urea. However, the degree of fluctuations 
in entropy and volume, characteristic of liquid H2O, is 
reduced [2j]. 

There was found [3(, that in water, two relaxation 
processes, differing only in activation entropy, may be 
due to two types of short-range order among molecules 
linked together by hydrogen bonds. Reorientation of the 
molecules occurs in structural defects, which allow the 
orientation of the molecules to change after the rupture 
a single hydrogen bond. The are two independent re- 
laxation processes correspond to that possible types of 
defects. 

In this work, we compare the behaviors of certain ther- 
modynamical values (an average energy and heat capac- 
ity) of aqueous AC and TMU in frames of statistical anal- 
ysis, which is useful for interpretation or confirmation of 
the hypotheses on mesomorphism [lj. In our modeling, 
we involve the Ising lattice version of the free energy of 
the systems, which was carefully studied in respect to nu- 
merical Monte Carlo methods [4|. In particular, binary 
water-surfactant systems are quite enough described 
[7[) via the low-temperature, the mean- field, and other 
approximations, and which are able to predict a major- 
ity mesomorphic states in aqueous solutions and can be 
tested numerically. 



ISING MODEL AND LATTICE MONTE CARLO 
SIMULATIONS 

The ability to formation of one or another space geom- 
etry of phases depends including on a geometry and size 
of surfactant molecules, i.e. on its volume fraction in a 
solution. Chemical properties of heads and tails of sur- 
factants, as well as their configuration thereupon, lead to 
formation of surfaces in solutions. The mesomorphism of 
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aqueous AC and TMU also obeys universal phase behav- 
ior in lyotropic transitions, described in 0, but a micellar 
degeneracy and other of these systems is outside of our 
report. 

The presence of some experimental values of chemical 
potentials of acetone at molar fraction 0.1 < x < 0.16, 
and for TMU at < x < 0.03 1], and in remainder 
ranges, calculated by the Boissonnas method, which al- 
low as to perform a calculations in the lattice model. 
Through there are the estimating formulas [H, 0] on re- 
lation the between fj, and x in aqueous AC, or a while 
below, we omit it for the checking of validity of a lattice 
modeling. 

Here, we apply the Ising gas lattice model Q , where an 
each lattice site u is assigned a spin variable a u = 0, 1, — 1 



of face-centered cubic lattice (fee). Depending on the 
molar fraction of particles, distributed on a such lattice, 
we have a vacations, labeled as —1, water molecules (0), 
and surfactant's molecules (AC and TMU), which have 
assigned a variables 1. We leave non-occupied lattice 
sites, because of arithmetic demand a total number of 
molecules in a system should be less or equal than a to- 
tal number of lattice sites, on which we distribute our 
system; so we rarely have their coincidence. An orienta- 
tional properties of surfactants are revealed via the coor- 
dinate variable s u and , where vectors s u are assigned 
each s ite of a lattice, and r vu specify the directions from 
fee lattice site u to v, for convenience, \r vu \ = 1. The 
Hamiltonian of the Ising model is: 



(1) 



r 



where < uv > denotes a summation over all distinct pairs 
of nearest neighbor sites u and v. The term with coupling 
Ji represents the isotropic interactions between nearest- 
neighbor sites. The terms with coupling constant J2 de- 
scribe the orientation-dependent interaction between a 
AC/TMU molecule and a neighboring water molecule. 
There is agreed that J2 > correspond to a surfactant 
orientation with its head pointing toward a neighboring 
water molecule. The last term gives the chemical poten- 
tial energy of the AC/TMU molecules. We remain here 
in the grand canonical ensemble, when the surfactant 
concentration is controlled by the surfactant (AC/TMU) 
chemical potential fi s , rather than in the canonical en- 



semble where the surfactant concentration is fixed. Fur- 
thermore bellow, we will use a normalized values of excess 
chemical potentials ^ e s x , instead of {i s -P>t x were calcu- 
lated and measured (as was reported in [lj and references 
therein) depending on the surfactant mole fraction, they 
are only free of the linear terms. Such a substitution 
will not affect on our average energy during Monte Carlo 
circles, since according to ([T]), the reciprocal (first and 
second terms) and 'concentration' (third term) spin vari- 
ables are summed up independently from each other. 

To simplify the calculations, the Hamiltonian (Ji} is 
normalized by the isotropic interactions J\ , and came to 
the form: 



J 



H = ~ ^ <J u (o u + l)a v (a v + 1) - y ^2 + 1 )( s « " r ™>) + i a v + !)( s u • r ra )j 
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Basing on the results of low-temperature expansion [7| , 
we used the calculated relations between j'2 and )T S at 
non-zero temperatures t (t — )cbT/Ji, where fc^ is the 
Bolzmann factor, and T is absolute temperature). In 
this case, the excess chemical potentials be considered 
coinciding with the standard ones. 



Fortunately, the data of and 0, cited in 111, allow 
us to employ the known values of /z s , connected with its 



molar fraction at T = 25°C as 



jUj = RT In — l -jr, 



(3) 



where Xi is the mole fraction of i-component in solution 
and p® is the vapour pressure of pure i, i means AC/TMU 
molecules, (see Table 1 and Table 2, squeezed from |l|). 

So we performed our simulations, assigning the num- 
ber of MC runs per a configuration 10 6 and the number 
of thermalization steps equaling 10 5 . To minimize an en- 
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TABLE I: Some values of a mole fraction of aqueous acetone 
and chemical potentials at 25°C [l). 

XAC ^AC 

1.000 0.000 

0.770 0.223 

0.580 0.580 

0.350 1.200 

0.223 2.300 

0.175 3.000 

0.100 4.000 

0.000 5.000 



TABLE II: Some values of a mole fraction of aqueous TMU 
and chemical potentials at 25°C [lj]. 

XTMU fJ-TMU 

0.900 0.000 
0.750 0.100 
0.600 0.300 
0.500 0.400 
0.320 0.500 
0.080 0.550 
0.070 0.650 
0.020 0.950 
0.015 1.200 

ergy of a configuration, we used the Metropolis algorithm 
building-in the Fortran program. 

As a trial result, to confirm the presence of a transi- 
tion states, corresponding to the Ising model [7j at high 
temperatures, we exhibit the temperature dependencies 
of specific heat (Figs. 1, 2) at j 2 — 1.5, which was given 
from the expansion of the free energy, corresponding to 
the fee conditions of the task. 

RESULTS AND DISCUSSIONS 

In this section, we report and discuss our numerical 
results based on the model expressed by (1), in order to 
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FIG. 1: The average heat capacity of aqueous acetone in de- 
pendence of temperature, obeying by MC simulations. 
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FIG. 2: The average heat capacity of aqueous TMU in de- 
pendence of temperature, obeying by MC simulations. 

illustrate mesomorphism of aqueous AC /TMU. 

Thermodynamics and structural schemes 

From the Gibbs energy equation, the extensive ther- 
modynamical quantities TS-pV remain after simulations, 
where TS is the partial molar entropy, V is the volume 
of a given binary system. At the constant volume, TS 
behaves pro rata a free energy. From Figs. 3, 4, respec- 
tive smoothness of the average energy, plotted for the 
aqueous AC and TMU systems, may be obvious. In the 
real systems, the transition from Mixing Scheme I to II 
starts at point X at x ~ 0.08 (for AC and TMU) and 
ends at point Y, x ~ 0.26 However, according to the 
computer simulations, for both systems, the specific heat 

< c > exhibits an absent of extremum between the points 
X and Y (Figs. 5, 6). In our MC tests, for aqueous AC, 

< c > has a sharp maximum in the vicinity of the 'II-IIF 
transition (at x ~ 0.93), and remains smooth in the re- 
gion of the 'I-IF transition. Fig. 7 shows the results of our 
MC tests, especially some typical disordered states at the 
Mixing Scheme II between the revealed transitions. Be- 
cause of irrelevant number of empty lattice sites, only two 
molecular variables and 1 are visible here. The picture 
is similar for aqueous TMU as well. Regions with clus- 
tered acetone and clustered water molecules are clearly 
seen. 

The correlation functions (Fig. 8 and 9 for the systems 
of AC— water and TMU— water respectively) reflect an 
order parameter, analogous to the magnetization in a 

usual Ising magnetic. So, for AC water, they show 

the temperature dependence taken at points of a molar 
fraction x corresponding to the mixing schemes. 

At low content of the solvent, 0.2, the order param- 
eter reveals a non-smooth deep phase transition, which 
promptly reaches a plateau asymptotic with appreciable 
fluctuations. At an increase of concentration of AC, ex- 
tremum indicated a phase transition is shifted to an in- 
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FIG. 3: MC results for the dependence of average energy on 
a molar fraction of aqueous acetone (AC), relative units. 



FIG. 4: MC results for the dependence of average energy (in 
relative units) on a molar fraction of aqueous TMU. 



crease of temperature, and an each curve becomes mono- 
tonic at a lower value of amplitude, than the foregoing 
one. 

For the ACwater system, our numerical results are 
comparable with such of [l2[ , noted in [l[ , for the concen- 
tration dependence of the excess partial molar volume in 
the 'I-IF transition point. 

For illustration, the results of the ACAC volumetric 
interaction function are reported in we note, that the 
presence of water in our system as additional futures of 
the diagram of Fig. 8. Analogously for aqueous TMU 
(Fig. 9), the correlation functions demonstrate three dif- 
ferent types of phase behavior, which correspond to three 
different mixing schemes. 



Experimental evidence of clusterization for aqueous 
axcetone 

For the 'I-IF transition in AC, the MC results have 
qualitative correspondence with the statements of Q , ob- 
tained from the dielectric relaxation measurements and 
chemical analysis, that this is just the phase transition, 
connected with the second relaxation process at x <0.4, 
at a reorientation of the water molecule, involving a dis- 
ruption of the inter-molecular hydrogen bonds in the 
dense structure. However, it is difficult to compare them 
directly because of different temperature, in which the 
chemical potentials were determined, and also thereupon 
that, in |l| and references therein, experimental excess 
chemical potentials were reported, while their total val- 
ues were used in [sj]. For the system 'water acetone', 
coincidence of the calculated above and reported in Q 
data for minimum of the free energy at x about 0.8 was 
revealed. 

Improvement neutron small-angle scattering measure- 
ments and theoretical treatments via the angle-averaged 
radial distribution functions technique Q allow us to 
compare the hypothesis by Koga with the real states of 
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FIG. 5: MC results for the dependence of average heat capac- 
ity (in relative units) on a molar fraction of aqueous acetone 
(AC). Error bar is not shown. 



molecular bonding in the case of the ACE water sys- 
tem. 

More exotic independent data on aqueous AC have 
been reached with positron annihilation spectroscopy 
13l 14 1, where it was confirmed not only the mixing re- 
gions of the Kogas studies, but the mechanism of clus- 
terization was proposed in frames of atomic interactions 
(hydrogen bonding and so on). 



CONCLUSION 

In polar aprotic solvents, the links with water are not 
strong, they are caused, in particular, by hydrogen bond- 
ing 3]. This paper answers the question about the role 
of this bonding in mesomorphism of the whole system. 
With regards to above statistical analysis, one may con- 
firm the data on dielectric relaxation [3| by revealing of 
the phase transition at high AC (and TMU) concentra- 
tions, and at once, received an absent of coincidences at 
the low-concentration region. Viz., at the lyotropic struc- 
ture rearrangement, there is not a distinct phase transi- 
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FIG. 6: MC results for the dependence of average heat capac- 
ity (in relative units) on a molar fraction of aqueous TMU. 
Error bar is not shown. 



tion in a system completely. Therefore, the expected 
structures stay beyond the Ising model description in the 
region of high molar fraction of AC and TMU, even if it 
possesses more sensitivity to a transition. In this prob- 
lem, it is important to separate a knowledge on pure 
water aggregates (hydrogen networks and cluster struc- 
tures [3, 4]) and on binary water mixtures (aqueous AC 
or TMU), comparison the behavior of hydrogen bond- 
ing in water with that in solutions may stimulate further 
advance in the theory of water. Other aqueous aprotic 
solvents might be also used to contribute such data. 
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FIG. 8: MC results for the temperature dependence of the 
correlation functions of aqueous acetone in relative units. 
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FIG. 9: MC results for the temperature dependence of the 
correlation functions of aqueous TMU in relative units. 



